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Abstract. We show that Casimir-force calculations for a finite number of non- 
overlapping obstacles can be mapped onto quantum-mechanical billiard-type problems 
which are characterized by the scattering of a fictitious point particle off the very 
same obstacles. With the help of a modified Krein trace formula the genuine/finite 
part of the Casimir energy is determined as the energy-weighted integral over the 
log-determinant of the multi-scattering matrix of the analog billiard problem. The 
formalism is self-regulating and inherently shows that the Casimir energy is governed 
by the infrared end of the multi-scattering phase shifts or spectrum of the fluctuating 
field. The calculation is exact and in principle applicable for any separation(s) between 
the obstacles. In practice, it is more suited for large- to medium-range separations. 
We report especially about the Casimir energy of a fluctuating massless scalar field 
between two spheres or a sphere and a plate under Dirichlet and Neumann boundary 
conditions. But the formalism can easily be extended to any number of spheres and/or 
planes in three or arbitrary dimensions, with a variety of boundary conditions or non- 
overlapping potcntials/non- ideal reflectors. 



PACS numbers: 03.65.Nk, 03.65.Sq, 03.70.+k, 05.45.Mt 
1. Introduction 

In 1948, the Dutch physicist H B G Casimir predicted the remarkable effect [1] that two 
parallel, very closely spaced, uncharged metallic plates attract each other in vacuum. 
The origin of this force can be traced back to the modifications of the zero-point 
fluctuations of the electromagnetic field by the presence of the two plates (at distance 
L) relative to the free case, or rather, relative to the case where the plates are separated 
by an infinite distance. Mathematically, this corresponds to the following difference of 
two mode sums: 

2 plates(L) ^ 2 plates(L^oo) 

The distinctive feature of the Casimir effect is that it depends on the geometry in a 
non-intuitive way: its strength and, perhaps, its sign are geometry-dependent (for a 
review see [2]). 

Let us look at the other side of the 'coin' and take the geometry dependence as 
guiding principle for the construction of the Casimir effect; i.e. we will define the 
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Casimir energy as that part of the vacuum energy that results from the geometry- 
dependent and therefore changeable part of the pertinent density of states. In fact, in 
Casimir-type calculations, the density of states p(E) = J^E k $(E — Ek), where [Ek] are 
the eigenenergies of the modes, can be split into three parts: 

p(E) = Po (E) + p bulk (£) + 6 Pc (E) . (1) 

(i) There is a free part Po(E) over the homogeneous part of the background (the 
fluctuating field, matter fields etc.) that is completely unaffected by the presence of 
the obstacles; (ii) the bulk part ptmik(-E) sums up the modifications of the spectrum 
(excluded volumes and surface contributions including Friedel oscillations) by each of 
the individual obstacles as though they were still infinitely separated; (iii) finally the 
remaining part Spc(E) is the genuinely geometry- dependent part of the density of states 
which 'knows' about the relative separations and angles between the various obstacles. 
It is the latter that determines the Casimir energy when it is integrated up, weighted 
by the energy E: 

roc rco 

Sc = \ dEE5 Pc {E) = -\ dEAf c (E) , (2) 

« " 

where Afc(E) = Jq dE' dpc(E') is the geometry-dependent part of the integrated density 
of states (or number of states), 

N{E) = [ E dE' p(E') = Y,Q(E - E k ) . (3) 

Note that the integrals in (T5]) do not converge absolutely, but only conditionally because 
of the oscillating behaviour of 5pc{E) and Mc{E). The inclusion of e.g. an exponential 
damping factor (compare with the analogous tilt of the integration path(s) in f[T3"j) ), 
which can be removed after the integration(s) are performed, cures this problem and 
also eliminates the spurious upper-boundary contribution ^EAf c (E)\°° from the partial 
integration. The lower boundary term \ENq{E)\q vanishes automatically. 



2. The mapping of the Casimir calculation onto a scattering problem 

The Casimir calculation for two parallel plates at distance L was simple since the 
problem is separable, i.e. quasi-one-dimensional if the plates had infinite extent. 
For more complicated geometries than two parallel plates only the proximity force 
approximation [3] is left in general (see, however, section |6]). 

However, as shown in [U [5] for the non-relativistic fermionic Casimir effect and in 
[HI [7] for the scalar Casimir effect^, the Casimir energy between any finite number of 
non-overlapping spherical obstacles, i. e. spheres and cylinders in 3 dimensions or disks 
in 2 dimensions, can be solved exactly, although the corresponding quantum-mechanical 
problems are not separable any longer. In fact, these calculations simplify because of 

| In the scalar Casimir effect the fluctuating field is not the electromagnetic one, but a massless scalar 
one with the dispersion E = hc\k\. In the non-relativistic fermionic Casimir effect, the background 
does not consist of fluctuating fields, but of (a Fermi sea of) non-relativistic fermionic matter waves. 
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Krein's trace formula [5] which maps - under the condition that the potentials of the 
obstacles are sufficiently short-ranged - the difference in the total level densities - with 
and without obstacles - to the energy variation of the total phase shiflQ. For the special 
case of n spherically symmetric obstacles of radii at mutual center-to-center separation 
Tij > a,i + dj (such that they do not overlap) and mutual angles a^- the Krein formula 
reads 

5p(E) = p{E) - p (E) = -L iL tr\n S n (E, {a,}, {%}) 

2m aE 

= 9^ ^Indet^^Ja,},^}), (4) 

where the distances and angles atij have been combined to separation vectors rij and 
where i,j = l,2,---,n label the obstacles. Note that the level densities on the left 
hand side are averaged over an energy-interval larger than the mean-level spacing in the 
volume V of the entire system which will be taken to infinity in the end. This averaging 
is done in order to match to the right hand side which is defined in terms of the total 
S'-matrix S n (E, {a«}, {rij}) of the n-disk system in 2 dimensions PJJ [TTJ [12] or the n- 
sphere system in 3 dimensions [13], respectively, and which is therefore continuous. In 
fact, in [Til [12] it was shown that the determinant of n-disk S-matrix is finite (except 
at the known poles and branch points), although the matrix has an infinite number 
of components, since the pertinent T-matrix was proved to be trace class. The 3- 
dimensional n-sphere case was discussed in [13]. In this way the Casimir calculation for 
a finite number of non-overlapping obstacles is mapped onto the quantum-mechanical 
analog of a classical billiard-type problem which is characterized by the scattering 
of a fictitious point particle off the very same obstacles. For the case of only one 
spherical obstacle the scattering problem is separable, but uninteresting as there is no 
Casimir effect. For two spherical obstacles the scattering problem is non-separable and 
classically hyperbolic, and for more than two obstacles it is in general even classically 

chaotic pn on na on in] . 

Let us remark that the energy-averaging and the infinite volume limit do not 
commute, but that first the volume has to be taken to infinity and only then the 
averaging interval can go to zero. This is obvious if one takes semiclassical considerations 
into account where the determinant over the n-sphere/disk matrix is given in terms of 
periodic orbits and possible Weyl term corrections, see [UJ [12]. For this purpose let 
us put the scattering system at the middle of a large container which can be taken, 
without loss of generality, to be spherical of finite radius R and consider this container 
together with its empty reference container that is of equal size. Now, as long as 
the containers have a finite size, there exist three classes of periodic orbits, (i) The 
generic scattering periodic orbits that only bounce between the obstacles belonging to 
the scattering system and that semiclassically specify the determinant of the scattering 
matrix, (ii) A first class of spurious periodic orbits which solely are there because of the 
container and which only bounce between the walls of the container. The subtraction of 

§ See |9] for a precursor of Krein's formula for potentials with spherical symmetry. 
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the reference system, in fact, removes this class of spurious orbits, (iii) A second class 
of spurious periodic orbits which bounce between the container walls and the scattering 
obstacles and which do not exist in the empty reference container. As the length L P . Q . 
of these periodic orbits will increase with the size of the container, their contribution 
to the periodic orbit sum can be suppressed if a small imaginary term ie is added to 
the energy or wave number k = \k\ which is here the same because the dispersion is 
E = hc\k\. The phase of the orbit acquires therefore an exponential damping factor 

which smooths (averages) the total and background densities of states in the Krein 
formula 

lini lim {p(k+ie, R) - p (k+ie, R)} = 1 . tr In S n (E(k), {ai}, {nj}). 

€-^0 + R-+00 ITU dhi 

This explains the smoothing/averaging procedure in (j3J). 



3. The mapping to the multi-scattering matrix 

In the following, we will show that the geometry-dependent Casimir fluctuations can 
be extracted from the mu/%>/e-scattering part of the S-matrix. For this purpose let us 
once more explain the philosophy: (i) we start out with a large, but finite container of 
radius R which will be filled with the background field (here the massless scalar field or 
a fermionic matter field) and which eventually will be taken to infinite size. The total 
density of states is therefore equal to the background density: p(k+ie, R) = p (k+ie, R). 
A suppression term is added to the wave number k as explained in the last section. 

(ii) Next we will also put the spherical obstacles of radii a« into the large container, 
where we make sure that the mutual separation r^ between any pair of obstacles is 
large, i.e. comparable with the size of the container. The total density of states has to 
be readjusted because the background field is affected by excluded volumes due to the 
obstacles and by boundary effects (including Friedel oscillations) due to the imposed 
boundary conditions at the surface of the obstacles, say Dirichlet boundary conditions: 

n 

p(k + ie, R, {ai}) = p (k + ie, R) + V p w (k + ie, R, a*) . (6) 

i=i' ' 

Weyl term 

(iii) Now, the obstacles are moved close to each other such that the actual system of 
a finite number of non-overlapping spherical obstacles is formed. Again the density 
of states has to be readjusted, this time by the geometry-dependent term 5pc which 
depends on all radii {a^} and on all separation vectors {r^}: 

n 

p(k + ie, R, {ai}, {%}) = p (k + ie, R) + Pw(k + ie, R, a*) 

i=l 

+ 5p c (k + ie,R, {a^}, {/%}) . (7) 

(iv) By inserting the difference of ([7]) and the background density of states p into the 
Krein formula and by taking the double limit lim e ^ 0+ h m R^oo we effectively remove the 



The Casimir effect as scattering problem 



5 



container walls and map the problem to the scattering system: 
lim lim {p((k + ie, R, {a*}, {/%}) - p {k + ie, R)\ 

= p(k, {a { }, {fij}) - p (k) = 7^-^ Indet S n (k, {a { }, {r^}) , (8) 

where lndet S n (k, {a^}, {?%}) is the total phase shift expressed through the n- 
sphere/disk S-matrix S n (k, {a,}, {?%})• As shown in [HI [12] and [13], respectively, the 
determinant of the n-disk/sphere S-matrix separates into a product of the determinants 
of the 1-disk/sphere S-matrices S\(E, ai), where czj is the radius of the single scatterer 
i, and the ratio of the determinant of the inverse multi-scattering matrix M and its 
complex conjugate: 

(» Wet (M{k*, {ai}, {fjj-})t) 

det S B (*, {a,}, {^}) = |n det S l (k, a,) J det M (fc, {aj, {f%}) " (9) 

Note that each and every determinant in ([9]) is finite and exists separately, since the 
pertinent T-matrices are all trace class [TTl [T2| [13]. 

When inserted into Krein's formula, the product over the single-scatterer 
determinants generates just the bulk (or Weyl term) contribution to the density of 
states 

n 

lim lim Pbuik(^ + ie, R, { a i}) = 5^ ^ m ^ m Pw(^ + ie, R, a i) 

e^0+ R— >oo i—Y e^0+ R^oo 

n 1 d n 

= z2pw(k,ai) = -— — 22ln det Sx(k,ai) , (10) 

i=l ^ i=l 

which takes care of the excluded volume terms and the surface terms (including Friedel 
oscillations). The geometry-dependent part of the density of states 

5p c (k, {ai}, {fij}) = lim lim 6p c {k + ie, R, {a*}, {r^}) (11) 

e — >U+ it — >oo 

is therefore given by a modified Krein equation [3] which is formulated in terms of the 
inverse multi-scattering matrix M {k, {ai}, {?%}) instead of the full S-matrix 

n 

Sp c (k, {a^, {r i:j }) = p(k, {a{}, [r i3 ■}) - p (k) - ^ p Wey i(&, a{) 

i=l 

= - -Im — In det M(k, {aj, {r%}) . (12) 

Note that in this expression all the possibly divergent terms connected with the sharp- 
surface limits are subtracted out. In fact, by mapping the system to a multi-scattering 
problem the calculation is self-regulating: an additional regulator is not needed (except 
the energy- damping mentioned below ([2]) which can be removed in the end, see also [6]). 

The pertinent Casimir energy can then be read off from the energy-weighted (or 
wave-number-weighted) integral 

/■oo /*oo 

S c = dE \E 5p c (E, {a^, {/%}) = \hc / dkk 5p c (k, {a*}, {f%}) 

J « 

he 



poo J^q roo 

\hc dk7V c {k, {<h},{rij}) = — / dk Im In det M(k) 

JO 2,71 Jo 
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he 

he 



oo(l+iO + ) 



dk In det M(k) 



oo(l-iO+) 



dfclndet M{k*y 



Tic r°° 

= — dfc 4 hidetM(i/c 4 ), (13) 

Z7T Jo 

where the dependences on the radii and separation vectors were suppressed in the 
argument of the inverse multi-scattering matrix M(k) = M(k, {a«}, {f^}) for simplicity. 
In the last step a Wick rotation k i/c 4 was performed in the first, and k — > — i/c 4 in the 
second integral, and the relation detM(i&4) = detM(i&4)' was applied which follows 
from the relation det M(k) = det M ((—&)*) [6]. As a corollary note that 



— I dk k 2n+1 lm\n det M(k) 



2tt Jo 



, he 

47X JO 



dki kj n + 1 



In det M{\k A ) - In det Af (ifc 4 ) 







(14) 



e.g. the Casimir energy over modes with a non-relativistic dispersion E = h /2mn, 
where mu is the pertinent mass, integrates to zero |6|, unless there is a finite upper 
cutoff, as e.g. the Fermi momentum fcp i n the case of the so-called fermionic Casimir 
effect 0]. 

Note that the final expression ffl3l) for the Casimir energy is obviously finite, since 
det M(ifc 4 ) vanishes rapidly with increasing value of fc 4 , and moreover, for the same 
reason, it is dominated by the infrared end of the integration. 



4. The multi-scattering matrix and the fermionic Casimir effect 

The modified Krein equation ([121) is especially useful as there there exists a close-form 
expression for the inverse multi-scattering matrix for n spheres (of radii aj and mutual 
distances r^y, where the indices j,j' = 1, 2, • • • , n label the spheres) in terms of spherical 
Bessel and Hankel functions of first kind, spherical harmonics and 3j-symbols [T3] : 

Mii Vm , = S jj '5 w S mm > + (1 - <P") i 2m+l '- 1 v /4tt(2Z+1)(2Z'+1) 

\a f ) hf{f\a f )korrk-v \ ) 

X ( m-m" I" -m )^U^h^(kr jf )Y^'{^). (15) 

Here 1,1' and m,m' are total angular momentum and associated magnetic quantum 
numbers, respectively. The rotation matrix D l m , m „(j,j') maps the local coordinate 
system of the sphere f to the one of the sphere j and the unit vectors f^} point from the 
origin of the sphere j (as measured in its local coordinate system) to the origin of the 
sphere j'. The above expression refers to Dirichlet boundary conditions on the spheres. 
The case of Neumann boundary conditions (or even mixed boundary conditions) on the 
spheres follows from the replacement 
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etc. A closed-form expression for the two-dimensional case of n disks, corresponding to 
(fl5j) . can be found in [IT] . 

For the case of small scatterers one can even simplify the expression of the multi- 
scattering matrix: 



e ifcr rf 



M» (k, {a,}, {r y }) « 5" - (1 - 5" )f s (k, a,) + 0(p-wave), (17) 

r jj' 

since the s-wave scattering dominates over all other partial waves. Thus, the non-trivial 
part of the inverse multi-scattering matrix is given by the propagation of spherical waves 
modulated by s-wave amplitudes f s (k,aj) between the spheres. From the modified 
Krein formula the integrated density of states in the case of two small spherical cavities 
of common radius a and center-to-center separation r can be deduced as |3] 

1 a 2 
Af°°(k) = ImlndetM°°(A;,a,r) » — sin[2(r - a)k] + O ((ka) 3 ) .(18) 

This expression should be compared with the semiclassical approximation that sums up 
all partial waves 

AfS%(k) = f_ sin[2(r - 2a)k}. (19) 



The latter is the leading contribution to Gutzwiller's trace formula [18] . namely the 
contribution of the non-repeated two-bounce periodic orbit between the two spheres. 
S po (k)/h = 2(r — 2a)k is the action of the two-bounce periodic orbit, where 2(r — 2a) 
is the length of its geometric path. Note that the semiclassical result is suppressed by 
a factor of 1/2 2 relative to the small-scatterer one. We will see that each factor 1/2 is 
associated with a semiclassical reflection from one sphere. 

It was shown in [3], under the condition k > 1/a, that the semiclassical result for 
Ac° sc (A;) is a very good approximation of the full quantum-mechanical result calculated 
from the exact expression ffT5]) of the two-sphere scattering matrix when plugged into the 
modified Krein formula (Tl2l) . Therefore, the Casimir energy for two spherical cavities 
inside a Fermi sea of non-relativistic non-interacting matter modes can be approximated 
in terms of a spherical Bessel function j\ as 

eg = -l k ' dm , a ,r),-^-^—^mr-2 a )k r] (20) 

where the expression is valid for kpa > 1. Here is the mass of the fermionic mode. 
Note that it is long-ranged, i.e. it scales as a 2 /L 3 with L = r — 2a. The corresponding 
fermionic Casimir energy of the sphere-plate system 

£c' - -^^7^ii[2(r - a)k F ] , (21) 
2m^s 2ir(r — a) 

even scales as a/L 2 with L = r — a. Moreover, the fermionic Casimir energy has in 
both cases an oscillating behaviour: with increasing distance between the obstacles, 
the Casimir energy, which starts out to be attractive, will become repulsive, and under 
a further increase of the distance, it will become attractive again, where the strength 
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is of course reduced. This is in contrast to the fixed (negative) sign of the standard 
Casimir effect with fluctuating electromagnetic or scalar fields between these obstacles. 
The reason for this oscillating pattern in the fermionic Casimir effect is the presence 
of a new scale, in addition to the length scale(s), namely the Fermi momentum /sp- 
in fact, the strength of this fermionic Casimir energy is determined by the UV-cutoff 
of the theory [I]. Also this behaviour distinguishes the fermionic Casimir effect from 
the standard Casimir effect: the latter is governed by the infrared behaviour of the 
corresponding density of states. 

In jl], the fermionic Casimir energy of the three- and four-sphere system was 
calculated as well. From the periodic orbit summation it is obvious that there exist 
genuine three- and more-body interactions. However, it was shown in [I] that the two- 
bounce orbit dominates in the equilateral three- and four-sphere systems. In fact, the 
billiard analogy holds: it is difficult to make long shots, especially with many bounces 
- the slightest error ruins the shot. Mathematically, the relative weights calculated 
from the instabilities of the longer periodic orbits, even weighted with their degeneracy 
factors, are far less than the weight of the two-bounce orbit. The maximal correction 
due to the next-to-leading periodic orbit, the triangular three-bounce orbit, is at most 
10 % at a distance r pa 2.5a and even smaller for other separations. 



5. The sphere-plate case of the scalar Casimir effect 

From now on we will only discuss the scalar Casimir effect, where the fluctuating modes 
are massless scalar ones. This effect was first investigated by Schaden and Spruch p^9| [20] 
for the case of two spheres subject to Dirichlet boundary conditions in the semiclassical 
framework of Gutzwiller's trace formula [18]. It was shown that the result of the 
proximity force approximation emerged from the semiclassical calculation in the limit 
of vanishing surface-to-surface separation. Gies and coworkers numerically studied the 
scalar sphere-plate case, again subject to Dirichlet boundary conditions, in the world- 
line approach [21]. Scardicchio and Jaffe analyzed these results in an optical approach 
that takes into account not only periodic, but also closed orbits [22J. 

In the case of two spheres there is a continuous axial symmetry with respect to 
the line joining the centers of the spheres, and an additional reflection symmetry with 
respect to any plane containing this (horizontal) symmetry axis. The symmetry group 
is Cooh in crystallography group theory notation [23], where the role of the vertical v 
plane and horizontal h plane is opposite to the one of this reference in order to match 
the orientation of figure [TJ As a consequence of the C^h symmetry the inverse multi- 
scattering matrix is separable with respect to the magnetic quantum number m: 

_ x ( ^lm,l'm I / 99 \ 

— °mm' I /l 21 x j \ I 

lm,l>m> \ A lm,l'm °W ) 

where the indices 1,2 label the two spheres and where M/^ z , m , is the left hand side 
of (fT5|) . whereas A^ mVrn i is the non-trivial term (minus the delta-distribution term 
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$i o~wo~mm') on the right-hand side of this equation - both adjusted to the two-sphere 
case. 

The case of two identical spheres at a center-to-center distance r contains the case 
of a sphere and a plate at a center-to-plate distance R = r/2 as can be seen in figured! 
Because of the additional two-fold reflection symmetry with respect to the vertical 



half-domain I 



vertical symmetry plane 



half-domain II 




Figure 1. Two spheres of common radius a at a center-to-center separation r. The 
vertical symmetry plane splits the full domain into two half-domains I and II. The 
center-to-plane separation R, the surface-to-surface separation L 00 in the full domain 
and L°l in the half-domain arc indicated as well. 



symmetry plane, the symmetry group of the two-sphere case with common radius a 
is even D^v- Therefore, the global domain splits into two half-domains (called I and 
II) that are separated by the vertical symmetry plane. Moreover, all scattering waves 
can be separated into symmetric ones (subject to Neumann boundary conditions) and 
antisymmetric ones (subject to Dirichlet boundary conditions) with respect to this plane. 
Thus there exist two classes of inverse multi-scattering matrices in one half-domain, one 
subject to Neumann boundary conditions (N) and one subject to Dirichlet boundary 
conditions (D) on the vertical symmetry plan [6]: 



Ma, 



oo (m) 



oo (m) 



(m) 



(using the short-hand notation A w 
determinants factorize as 



(m) 



A 12 

Im.l'm 



A 21 \ 



(23) 

and the associated 



det M°°(k, a, r) = f[ det M oo(m) (A;, a, r) 

m=— oo 

= det M°°{k, a, r)\ N det M°°{k,a,r)\ D . (24) 

Therefore the determinant for the sphere-plate case under Dirichlet boundary conditions 
can be inferred from det M°°(k, a, r)| D as follows: 

det M o1 (k, a, L g| ) = det M°° (k, a, r = 2(L g| + a)) | (25) 

and the pertinent Casimir energy is given as |6] 



£$(a,L) = — / dk 4 \ndet M°°(ih,a,r = 2(L + ( 
2n Jo v 



D 



(26) 
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The exact data generated by this equation are reported in [5] . They are compatible 
with the numerical data of the worldline approach [21] between L = a/8 and L = a if the 
statistical error bars of the worldline data are taken into account, see [6]. Moreover, the 
(systematically and statistically improved) worldline data of [231 E5j , which extend up 
to L = 16a and up to L < 100a, respectively, do nicely agree with the exact data of [6]. 
In principle, the calculation based on (126]) is exact and applicable for all separations of 
the obstacles. In practice, it is more suited for large- to medium-range separations. The 
reason is that the Casimir energy for the scalar case is dominated by momenta k ~ 1/L 
where L is the separation scale. As shown in [6] for the scalar sphere-plate case, the 
Casimir integration can be truncated at k max ~ 10/L corresponding to a truncation in 
the angular momentum space at Z max < (e/2)fc max a « 14a/ L. For L < 0.1a the matrices 
increase rapidly in size and the angular- momentum algebra (see the 3j symbols in ( fl5l) ) 
becomes very cumbersome. 

Let us rather concentrate on the asymptotical behaviour (surface-to-surface 
separation L ^> a) of Casimir energy given by ( 1261) . The s-wave approximation is 
not needed for the fermionic Casimir energy since it is governed by the UV part of 
the density of states (i.e. by the contribution at the Fermi momentum fcp, assuming 
k-p > I /a). However, it dominates the asymptotics of the ffuctuating-scalar Casimir 
effect at very large separations (L ^> a), as this Casimir-energy type is governed by the 
infrared behaviour of the density of states at k ~ 1/L or less. The Casimir energy of the 
sphere-plate case subject to Dirichlet boundary conditions is given at large separations 
L ^> a (in the s-wave approximation) as [6] 

7T 3 hca 90 2 

^ ' ~ ~ 1440L 2 7^(1 + a/L)(l + a/2L) 

7r 3 hca 90 n 3 hca , _ s 

- -l4i0Z^ X2 = -l440F X 1 - 847 "" < 27 > 
As in section HJthere is a relative factor of 2 originally found in [3] between the s-wave and 
semiclassical result for each reflection off a sphere. This applies also for the fluctuating- 
scalar Casimir effect, if the relativistic dispersion E = Tick and the suppression of the 
repeats of the semiclassical two-bounce orbit at large separations (by the removal of 
the term Xi n =i w" 4 = 7r 4 /90) are taken into account [6]. Therefore, the exact result of 
the fluctuating-scalar Casimir energy for the very far separated Dirichlet sphere-plate 
system is enhanced by a factor of 2 relative to the semiclassical result and by a factor 
of 2 x (90/7T 4 ) relative to the leading term — 7r 3 aca/(1440L 2 ) of the proximity-force 
approximation (PEA) [3]. The leading term of the PFA, in other words the PFA at 
vanishing distances, where the repeats of the two-bounce orbit cannot be neglected and 
add up to the factor 7r 4 /90, was - as mentioned - confirmed semiclassically in [19J. In 
this connection, it should be noted that the result quoted in [6J for the semiclassical 
approximation to scalar Dirichlet sphere-plate case 



he a ( 7T 



16tt L 2 \ 90 



5 1\ L 



(28) 
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also applies for the scalar Neumann case, since both semiclassical calculations differ only 
by the Maslov indices ii.e by a minus sign for each reflection from a Dirichlet surface) 
and since the number of reflections is even for a self- retracing orbit. The electromagnetic 
case in the semiclassical approximation is then given by the sum of both scalar cases. 
It is therefore twice as big as (128j) and the first correction to the leading PEA result is 
predicted to be negative. 

In the case of the electromagnetic Casimir effect for the sphere-plate system, the s- 
wave dominance at large separation has to be replaced by a p-wave dominance, since the 
charge- neutrality of the sphere forbids a monopole term, whereas the standard Casimir- 
Polder energy is dominated by induced-dipole contributions. If one removes by hand the 
s-wave contribution from the Casimir energy of the scalar Dirichlet sphere-plate system, 
the remaining energy is dominated at large separations by the p-wave contribution 

57r 3 hca 3 90 

£>wavc(£)~- 1440L 4 ( 29 ) 

Note that this expression is compatible with the a 3 /L A scaling of the Casimir- Polder 
energy for a molecule-plate system [26], but the prefactor is different [27]. 

Together with M Bordag [28] the asymptotics of the Dirichlet sphere-plate problem 
was worked out to sixth subleading order in the ratio of the sphere radius a and the 
center-to-plate separation R — r/2: 

hca r 5 a 421 / a \ 2 535 (a \ 3 3083041 



a 



8ttR 2 { 8R lUKRj 1152 V-Ry 518400 \R 

2741117 / a \ 5 557222415727 (a \ 6 1 , , 

- +••• . (30) 



1382400 \RJ 36578304000 \R, 
Even higher-order contributions (up to ninth order) were considered in this asymptotic 
series, but did not lead to any improvements. For practical purposes the sixth order 
expression is still useful down to separations L = 0.7a (or R = 1.7a), see figure [21 If the 
s-wave contribution is removed by hand, the corresponding expansion reads [28] 

5a 3 f 56 fa\ 2 597 /a\ 3 10453 /a\ 4 



16nR A I 25 \RJ 640 \RJ 1750 \RJ 
16557 / a \ 5 394844679647 fa \ 6 1 

o + M . CTe nnn n +"" " ( 31 ) 



1600 \RJ 9144576000 \R, 
Finally, for the corresponding Neumann sphere-plate system the removal of the s-wave 
contribution is natural and the corresponding result is given by [28J 

10a 3 f 63 fa \ 2 597 / a \ 3 4159 /a\ 4 



En,1>0 167T.R 4 L 1 + 100 \RJ + 320 \RJ 14000 \RJ 
271437 fa \ 5 148355331834 / a \ 6 1 

25600 \RJ + 2286144000 \RJ + " 7 ' ^ ' 

For reasons of completeness, we also report the numerical result of the pure s-wave 
contribution for the sphere-plate case under Neumann boundary conditions, although 
it is of course artificial: 

E^i=o = x 0.46066 . . . x jl + O (-|) } . (33) 
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Figure 2. Exact and asymptotic results for the scalar Casimir energy of the sphere- 
plate configuration with Dirichlet boundary conditions are shown shown in units of 
— fi.c7T 3 a/(1440£ 2 ) as function of the ratio of the surface-to-surface-separation L and 
the sphere radius a. The exact result from (|26[) is plotted as full line (in red), the 
asymptotic approximation (|30[) up to the (a/R) term as dashed line (in green) and 
up to the (a/R) 6 term as short-dashed line (in blue). 



Note that in this case the infrared limit — > and the asymptotic limit R/a — > oo do 
not commute in ( J26l) . 

Finally, for comparison, the result of the Casimir-Polder calculation in the limit of 
a perfectly conducting sphere as calculated in [27] is listed here as well: 



We see that the relative ratio of the leading Dirichlet and Neumann I > contributions 
is correct, but that the scalar result is bigger than the electromagnetic one by a factor 



6. Final remarks 

In the meantime, after the work [6] appeared, further exact Casimir results for non- 
separable systems have been reported. In [29], the Casimir interaction between a plate 
and a cylinder was calculated for the background of a fluctuating electromagnetic field. 
The exact zero-point interaction between two non-concentric cylinders was studied in 
[30] . The Casimir effect for a sphere and a cylinder in front of a plane and corrections 
to the proximity force theorem were investigated in [31J. The Casimir forces between 
arbitrary compact objects were investigated in |32j; see also the study [33] in the T 
operator approach and pMj [35]. 




(34) 



of 5/3. 
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Finally, let me mention the way in which the Casimir energy per unit length for 
n non-overlapping parallel cylinders of infinite length can be derived from the modified 
Krein formula of the two-dimensional n-disk case. For a cut perpendicular to the 
cylinder axes the modified Krein formula of the n-disk case can be applied, provided the 
wave number k in that formula is replaced by the modulus of the perpendicular wave 
number k± (the wave number along the cut). Moreover, the phase space integration 
L|| / d/c||/(27r) along the parallel direction to the cylinder axes has to be added and the 
resulting integrals have to be weighted by the correct energy-dispersion hc^Jk^ + k\. 
The corresponding Casimir energy per unit length therefore reads 

/OO dk\\ /"OO , 
— — / dk\_\%cJk\ + k\8pc{k±, n-disk) 
-oo 2n Jo v 11 

/oo dk\\ roo qI — \ 
— - / dk ± ±hcjk? l +k 2 l — ImlndetM(/c_L,n-disk) 
-oo 2n Jo * 11 dk± 7r 

= — 7T ( d/cii / dk 4 lndet M(i/c 4 , n-disk) 

2W-00 ll J\ kll \ \JkJ^q 1 ; 

[•OO 

I dk 4 k 4 In det M(i/c 4 , n-disk) . (35) 
Jo 



he 

An 



7. Summary 

We have shown that the Casimir energy can be re-defined as the vacuum energy 
of the geometry-dependent part of the density of states and that the latter can be 
calculated form the multi-scattering phase shift of a modified Krein formula. The non- 
overlapping, in general non-separable n-sphere, sphere-plate, n-disk (and n-cylinder) 
Casimir problems can be solved exactly in the scalar (and also in the fermionic) case. 
The calculation is not plagued by the subtraction of single-sphere contributions or by 
the removal of diverging ultraviolet contributions. All involved determinants exist and 
are finite since the pertinent T-matrices are trace-class. The large- distance behaviour is 
dominated by the s-wave scattering in the case of the scalar Casimir effect with Dirichlet 
boundary conditions and by the p-scattering for the corresponding Neumann case as in 
the electromagnetic scenario. The Dirichlet or Neumann boundary conditions can even 
be replaced by mixed boundary conditions. The presented method can easily be applied 
to any number of spheres or cylinders with or without planes (in two dimensions: disks 
with or without lines). Moreover, the spheres (or disks) can be replaced by other objects 
or even smooth potentials or non-ideal reflector, as long as these objects do not overlap. 
The finite surface thicknesses can be booked as Weyl-term contributions. 



Acknowledgments 



A.W. would like to thank the organizer Michael Bordag for the invitation to this excellent 
workshop QFEXT'07 in Leipzig. 



The Casimir effect as scattering problem 



14 



References 

Casimir H B G 1948 Proc. Kon. Ned. Akad. Wetensch. 51 793 
Bordag M, Mohideen U and Mostepanenko V M 2001 Phys. Rept. 353 1 

Derjaguin B V, Abrikosova I I and Lifshitz E M 1956 Quart. Rev. 10 295; Blocki J, Randrup J, 

Swiateck W J and Tsang C F 1977 Ann. Phys., NY 105 427 
Bulgac A and Wirzba A 2001 Phys. Rev. Lett. 87 120404 |arXiv:nucl-th/0102018j 



Bulgac A, Magierski P and Wirzba A 2005 Europhys. Lett. 72 327 arXiv:cond-mat/0406255 



Bulgac A, Magierski P and Wirzba A 2006 Phys. Rev. D 73 025007 farXiv:hep-th/05 1 1056| 
Wirzba A, Bulgac A and Magierski P 2006 J. Phys. A: Math. Gen. 39 6815 
arXiv:quant-ph/0511057| 



Krein M G 1953 Mat. Shorn. (N.S.) 33 597; Krein M G 1962 Sov. Math.-Dokl. 3 707; Birman M 

Sh and Krein M G 1962 Sov. Math.-Dokl. 3 740 
Beth E and Uhlenbeck G E 1937 Physica 4 915; Huang K 1987 Statistical Mechanics, John Wiley 

k Sons, New York, ch. 10.3; Friedel J 1958 Nuovo Cim. Ser. 10 Suppl. 7 287 
Gaspard P and Rice S A 1989 J. Chem. Phys. 90 2225; 90 2242; 90 2255 
Wirzba A 1999 Phys. Rept. 309 1 |arXiv:chao-dyn/9712015| 

Wirzba A and Henseler M 1998 J. Phys. A: Math. Gen. 31 2155 [arXiv:chao-dyn/9702004 



Henseler M, Wirzba A and Guhr T 1997 Ann. Phys., NY 258 286 |arXiv:chao-dyn /9701018 

Eckhardt B 1987 J. Phys. A: Math. Gen. 20 5971 

Cvitanovic P and Eckhardt B 1989 Phys. Rev. Lett. 63 823 

Wirzba A 1992 CHAOS 2 77 

Wirzba A 1993 Nucl. Phys. A 560 136 

Gutzwiller M C 1990 Chaos in Classical and Quantum Mechanics (Springer, New York) 
Schaden M and Spruch L 1998 Phys. Rev. A 58 935 
Schaden M and Spruch L 2000 Phys. Rev. Lett. 84 459 

Gies H, Langfeld K and Moyaerts L 2003 J. High Energy Phys. 06 018 |arXiv:hep-th/0303264] 
Jaffe R L and Scardicchio A 2004 Phys. Rev. Lett. 92 070402; Scardicchio A and Jaffe R L 2005 
Nucl. Phys. B 704 552 

Hamermesh M 1962, Group Theory and Its Applications to Physical Problems (Addison- Wesley, 
Reading, MA), Chap. 9-1 



Gies H and Klingmuller K 2006 J. Phys. A: Math. Gen. 39 6415 |arXiv:hep-th/0511092| 



Gies H and Klingmuller K 2006 Phys. Rev. D 74 045002 |arXiv:qua nt-ph/0605141 
Casimir H B G and Polder D 1948 Phys. Rev. 73 360 
Datta T and Ford L H 1981 Phys. Lett. 83A 314 
Bordag M and Wirzba A unpublished 

Emig T, Jaffe R L, Kardar M and Scardicchio A 2006 Phys. Rev. Lett. 96 080403 

arXiv:cond-mat/0601055| 
Mazzitelli F D, Dalvit D A R and Lombardo F C 2006 New J. of Phys. 8 240 

|arXiv:quant-ph/0610181] 
Bordag M 2006 Phys. Rev. D 73 125018 |arXiv:hep-th/0602295| 

Emig T, Graham N, Jaffe R L and Kardar M 2007 larXiv:0707. 1862 [cond-mat.stat-mech] and 

larXiv:0710.3084l [cond-mat.stat-mech] 
Kenneth O and Klich I 2007 larXiv:0707.4017l [quant-ph] 
Emig T and Jaffe R L 2007 larXiv:0710.5104l [quant-ph] 
Milton K A and Wagner J 2007 larXiv:0711.0774l [hep-th] 



